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Abstract 

We relate the relativistic finite range mean-field model (RMF-FR) to the point- 
coupling variant and compare the nonlinear density dependence. From this, the 
effective Hamiltonian of the nonlinear point-coupling model in the nonrelativistic 
limit is derived. Different from the nonrelativistic models, the nonlinearity in the 
relativistic models automatically yields contributions in the form of a weak density 
dependence not only in the central potential but also in the spin-orbit potential. 
The central potential affects the bulk and surface properties while the spin-orbit 
potential is crucial for the shell structure of finite nuclei. A modification in the 
Skyrme-Hartree-Fock model with a density-dependent spin-orbit potential inspired 
by the point-coupling model is suggested. 

Key words: Skyrme-Hartree-Fock model, relativistic mean-field model, 
nonrelativistic limit 
PACS: 21.30.Fe, 21.60.Jz 



1 Introduction 

Relativistic point-coupling (RMF-PC) models have proven to deliver predic- 
tions for nuclear ground-state observables which are of comparable quality as 
the ones from the well-established finite-range relativistic mean-field (RMF- 
FR) model (for a review see [4]) and the Skyrme-Hartree-Fock (SHF) approach 
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(for a review see [6]) [1,2]. Besides opening up the way to relativistic Hartree- 
Fock calculations (see Ref. [7] for a recent application) numerically similar to 
Hartree calculations with the use of Fierz relations for the exchange terms 
(up to fourth order, see Ref. [3]) and to the study of the role of naturalness 
[8,9] in effective theories for nuclear structure related problems, the RMF-PC 
approach also provides an opportunity to study the interrelations between 
nonrelativistic and relativistic point-coupling models, i.e., between the RMF- 
PC and the SHF approach. In Ref. [22] Rusnak and Furnstahl have shown 
the profitability to apply the concepts of effective field theory such as natural- 
ness to point-coupling models, where besides the version which will be used in 
our analysis, they consider also tensor terms, the mixing among the densities 
in the nonlinear terms, and nonlinear derivative terms. This investigation is 
motivated by the fact that up till now the role and the importance of the var- 
ious terms in the RMF-PC ansatz is not completely understood, not speaking 
about terms which might yet be missing. There appear systematic differences 
in the model predictions which could not yet be mapped onto the correspond- 
ing features of the models. A comparison between the nonrelativistic limit of 
the RMF-PC model and the SHF model may help to clarify these questions. 

One problem to face when comparing relativistic finite range with point- 
coupling nonrelativistic models is that two limits have to be taken: (a) the 
limit of letting the range of the mesons shrink to zero and (b) the expan- 
sion in powers of v/c (nonrelativistic reduction). The connection between 
the Skyrme-Hartree-Fock model and the RMF-FR model was done by sev- 
eral authors, without [10,11] and with nonlinear terms [4,12] (employed as 
self- interactions of the a-meson), but they did not take into account tensor 
contributions. The role of the tensor coupling of the isoscalar vector meson 
to the nucleon in the framework of effective field theories was investigated in 
Ref. [13]. 

The nonlinear density dependence is introduced in much different fashion for 
Skyrme-Hartree-Fock as compared to the Walecka model which employs non- 
linear meson self-couplings. This is a hindrance for a direct comparison [4,12]. 
On the other hand, nonlinear terms in RMF models are important, because 
only relativistic models with nonlinear terms can reproduce experimental data 
with acceptable accuracy [4,14,1]. We avoid the problems if we use the RMF- 
PC model [1], because in this model the nonlinear terms are explicitly density 
dependent, similar as in Skyrme-Hartree-Fock. Therefore it is worthwhile to 
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study the connection between RMF-FR and RMF-PC on one hand and be- 
tween the RMF-PC model and the nonrelativistic Skyrme-Hartree-Fock model 
on the other hand. In this context, we can study the role not only of the linear 
terms but also the nonlinear ones of both models in the nonrelativistic limit. 

The paper is outlined as follows: Section 2 evaluates the expansion of the 
finite range meson propagators of the RMF-FR into point-couplings. The non- 
relativistic limit is then discussed in section 3 from the RMF-PC as starting 
point. The part discussing in particular the emerging structure of the spin- 
orbit functional is taken up in section 4. And a few general comments on 
exchange are finally made in section 5. 



2 From RMF-FR to RMF-PC 



The RMF-PC model can be considered as the mediator between RMF-FR and 
SHF. The effects of finite range have nothing to do with the nonrelativistic 
limit, so we study them first by comparing the zero-range limit of RMF-FR 
with the point-coupling ansatz while remaining at the level of RMF. Having 
done this, we proceed in the subsequent section with the derivation of the 
nonrelativistic limit of the RMF-PC model. 

The covariant formulation of the RMF is based on a Lagrangian density. It 
is given for the RMF-FR in appendix A. 2 . For the stationary case, it can 
equally well be formulated as a Hamiltonian density. This reads 

7i = Tifrcc + T^s + + Wr, (la) 
C = E *«H7- V + m B )* Q , (lb) 

a 

Ks = \ ((V$) 2 + m|$ 2 ) + gs&ps + \b 2 e + i& 3 $ \ (lc) 
H Y = ~ (VPVV, + m 2 v V^) + gypoVo - ^PrV , (Id) 
n R = - l - (VR^VR^ + m 2 v R»R^ T ) + g R p T0 R T0 - \ ^P^RtO ■ 

(le) 

<7i, fi are coupling constants and the indices i denote scalar (S), vector (V), 
tensor (T) and isovector (r). $, V and R r0 are the isoscalar-scalar and the 
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zero components of the isoscalar-vector and isovector-vector meson fields, re- 
spectively. The densities are defined as corresponding local densities: 



isoscalar-scalar: ps(j) = J2 a 4>a if)^ (^), 

isoscalar-vector: p (r) = £ Q 0q(^7o0q(O> 

isovector-vector: p T o(r) = E Q 4>a(r)T 3 'yo<f>a(r) , ( 2 ) 

isoscalar-tensor: p^if) = — iE«V- ((p a (r)ci(f>a{r)), 

isovector-tensor: p T T{r) = — iE«V- (0 Q ,(r)r 3 a0 Q ,(f)). 

Now the Hamilton density will be expressed exclusively in terms of the densi- 
ties (2). To this end, the meson fields are eliminated by inserting the solution 
of the meson field equation. This is straightforward for linear coupling. We 
exemplify it here for the isoscalar-vector field. The meson field equation is 



(m v - A) Vq = g v p - t^Pt ■ 
This is solved for V by expansion in orders of A n going up to first order: 
9vV = 



1 - A/m 



9v 9vfv 1 _ 9v \ 

~ a v p - — PtH r A Po 

2 m v 

O-V 6 _T_ <V 



v_ 

rriy ru IrriBrriy 



2 Po ~ —_ .__ 2 Pt 
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Reinserting that into the Hamiltonian densities (Id) yields 

7-L v = ^a v pl + ^SyplApo - yPrPo + \^tP% (3a) 
with recoupled strengths 



a v = (3b) 
m v 

<V = 4" ( 3c ) 
m v 

a 9vfv 

9 T = J" ( 3d ) 

fy /o \ 

The form of the recoupled strengths is similar for the (linear) isovector-vector 
term. 
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The expansion is more complicated for the nonlinear isoscalar-scalar term. It 
involves a combination of expansion and iteration: first, the zero-range expan- 
sion as above, and second, an iteration of the nonlinearity. We start from the 
equation determining the scalar field 



9s® = 1 _ A ^ m2 (asps + b 2 ® 2 + h® 3 ) , (4) 

where as is equal to —gs/ m s an d 6fc is equal to —bugs I wis- The meson prop- 
agator is expanded to ~ 1 + A/m| as in the vector case. The nonlinearity is 
resolved by an iteration process. We can also obtain the $ by using another 
way, for example see Appendix A. 3. We obtain the structure 

Ks = ^a s pl + ^5 s p s Aps + ^(3sp s + \lsp s + CsPs^ps + C s p s Aps + .... (5a) 
with coefficients 



Ps = -b2^ , (5b) 

75 = (684-262-4) » ( 5c ) 
rrr s m s u 

Cs = -3& 2 4- > ( 5d ) 

<s=(^A-^A) ■ (5e) 



At this point, we can discuss the formal structure of the emerging effective 
Hamiltonian in comparison with the point-coupling model (6). We see that 
all the terms of RMF-PC are nicely generated by the above expansion. These 
are the terms with the coefficients a m , (3s, 75, S m , and 6 rn . The expansion 
generates some more terms not contained in the RMF-PC model. There is 
the term in K m , the finite range correction for the tensor term. It can be 
assumed to be small because the tensor coupling as such is already a small 
correction. And there are the many further terms generated by the expansion 
of the nonlinear $ coupling plus finite-range corrections thereof, i.e. the terms 
in £g, (' s etc. They require a more quantitative consideration. 

The mapping of the coefficients for known parametrisations can be seen in 
table 1. The first three columns show the effective point-coupling parameters 
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from the RMF-FR forces NL1 [4], NL3 [5] and NL-Z2 [21]. They are compared 
with the parameters of two genuine point-coupling forces PC-LA [20] and 
PC-F1 [2]. The "simple" parameters a m and 5 m agree nicely amongst the 



Parameter 


NL3 


NL1 


NL-Z2 


PC-F1 


PC-LA 


as 


-15.74 


-16.52 


-16.45 


-14.94 


-17.55 




-2.37 


-2.65 


-2.65 


-0.63 


-0.64 


OLy 


10.53 


10.85 


10.66 


10.10 


13.34 


S V 


0.67 


0.67 


0.68 


-0.18 


-0.17 


^TJ 














0.029 



















ce T V 


1.34 


1.65 


1.39 


1.35 


1.27 




0.09 


0.11 


0.11 


-0.06 





Ps 


38.11 


52.64 


58.74 


22.99 


3.32 


IS 


-347 


-590.77 


-706.63 


-66.76 


131.80 


Cs 


17.23 


25.37 


28.21 








Cs 


-196.78 


-348.92 


-408.00 



















-8.92 


-100.84 



Table 1 

Comparison of the point-coupling parameters between RMF-FR and RMF-PC 

various models. Comparing the S m between the two models, we notice that 
their absolute value is smaller in the point-coupling approach. Furthermore, as 
was already pointed out in Ref. [2], the 5 m values for the vector channels (both 
isoscalar and isovector) have a different sign than the ones from the RMF-FR 
variant. This strongly indicates that their role goes beyond the expansion of 
the propagators. 

The parameters associated with nonlinearity show large deviations. Moreover, 
one sees that the expansion of the nonlinearity within the RMF-FR approach 
is slowly converging. The expansion was done here around p$ = 0. One may 
hope that other expansion points, as e.g. bulk equilibrium density, lead to 
better convergence. We have checked that and find that the situation remains 
as bad. A detailed explaination can be found in appendix A. 3. A quick check 
is to take a typical value for the scalar density in bulk, p s ~ 0.14 fm~ 3 and 
to multiply each term with its power of ps- We have then a sequence of 0.32, 
0.14, 0.22 for NL1 and similar for NL-Z2. The terms in Cs are n °l small ei- 
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ther. This demonstrates that the parametrization of nonlinearity is different 
in its structure in both models. On the other hand, NL-Z2 and PC-F1 pro- 
duce very similar results for a broad range of observables in existing nuclei [2]. 
Actual observables explore only a small range of densities around bulk equi- 
librium density, and they do not suffice to assess the underlying differences in 
nonlinearity. 

To summarize: the expansion of the meson propagator of the RMF-FR into 
derivative couplings in RMF-PC works fairly well for the leading parameters. 
The resulting derivative couplings from RMF-FR are different in strength and 
sign from those of the RMF-PC which hints that there are other, more genuine, 
sources for gradient terms (quite similar to the density functional theory for 
electrons [26]). The worst case is the expansion of the non-linearity. The forms 
are so different in RMF-FR and RMF-PC that we could not find a simple 
mapping. In order to make these differences visible in practical applications, 
one needs yet to look for observables which are sensitive to very low densities. 
Halo nuclei could be a promising tool in that respect [27] . 

3 From RMF-PC to SHF 

In this section we study the nonrelativistic reduction starting from the RMF- 
PC model including tensor terms and nonlinear terms in both isoscalar-scalar 
and isoscalar-vector densities. As starting point we use the energy density of 
the RMF-PC model 



This Hamiltonian contains, besides the tensor terms (isoscalar and isovector), 
the isovector-vector term which appeared to be the most important one in 
former investigations [1,2]. The parameters 7$, 6i are usually determined 

in a x 2 adjustment to finite nuclear observables. The tensor terms can either 
be put in by hand in a Hartree theory (like in [22]) or be thought of emerging 




(6b) 



(6a) 
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from an approximate treatment of derivative exchange terms [7]. 

We now want to derive the nonrelativistic limit of that energy density following 
the procedure as described in [4]. In the first round, we consider only isoscalar 
fields and drop the Coulomb interaction to keep notations simple. Isovector 
contributions will be taken into account later in the calculation of the spin- 
orbit terms, because only in this sector the effect is significant. The isoscalar 
part of this Hamiltonian density leads to the stationary Dirac-equation 

[-17 • V + m B + S + 70^0 + id ■ f]* Q = e a7o ^ a , (7a) 

where 



S = a s ps + S s Ap 5 + (3sP 2 s + IsPs » ( 7b ) 

O'Y' 

V = a v p + 5 V Ap - —p T , (7c) 

T=-fV Po • (7d) 

The details of the nonrelativistic expansion are given in appendix A. 4. The 
nonrelativistic expansion in orders v jc oc p/m up to (p/m) 2 requires, of course, 
small p/m. One also needs to assume e a ~ which, however, is related to 
the first assumption of small momenta. The result of the expansion is that the 
normal nuclear density is associated with the zeroth component of the vector 
density p . We use henceforth the identification p = P- The scalar density 
and the vector density are expressed in terms of this nuclear density p and 
further densities and currents as follows: 



p s = p-2B 2 (r-V-J + pf 2 -2f-J + f-fp 

p T = -V • (B Vp) + 2V • (B J) - 2V • (B oP f), 
B = [2m B + S-V o r 1 . 

The nonrelativistic densities p, r and J are defined as 



(8a) 

(8b) 
(8c) 



a 

^ = E^(Vy? clt )-(Vv? cl ) , 



^ clt (v x oV cl ) - (v x av cl )V 



;;,„cl\t cl 



(8d) 
(8e) 

(8f) 
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where <p is the nonrelativistic single-nucleon wavefunction. 

Finally, we insert the expanded ps and pr into the energy density (6) keeping 
again terms only up to second order. This yields the nonrelativistically mapped 
energy density as 

^ (cl) = \{<*s + «vV + ^Psp 3 + \lsp 4 + \{5s + 5 v )p^P 
- (asp + Psp 2 + lsp 3 ) T 

- yP ("V • (B Vp) + 2V • (S J) + T V ■ (B oP Vp)) , (9) 

where 

T = 2B 2 Q (t - V • J + pf 2 -2T-J + T- Vp) 

= 2S 2 (r - V • J + ^p(Vp) 2 + ^Vp • J - y(Vp) 2 ) . (10) 

For better comparison, the Hamiltonian density is cast into a general form 

where basic structures are singled out with a separate coefficient each and the 
coefficients all may depend on the density, i.e. Cj = Cj(p). Less simple forms 
are lumped together in 5TL. This form can be directly compared with the 
standard nonrelativistic mean field model, the Skyrme-Hartree-Fock (SHF) 
energy functional, for a recent review see [24]. The energy functional of SHF 
is given in appendix A.l. It is obvious that it has the structure of the functional 
(11). Table 2 compares the coefficients for the here derived nonrelativistic limit 
of RMF-PC with SHF. The table shows the similarities and the differences 
between the two models. The striking similarity consists in the fact that all Cj 
terms appear in both models. A difference appears in the relation between r- 
and VJ-term. The RMF-PC without tensor coupling predicts C3 = C4 while 
SHF has two separate (and practically different) coefficients for that. It is 
interesting to note that the tensor coupling in RMF-PC also allows separate 
adjustment of these two terms. A basic difference appears with respect to the 
density dependence of the coefficients. All coefficients of RMF-PC carry a more 
or less involved density dependence. Only C\ is density dependent for SHF and 
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RMF-PC basic RMF-PC tensor 


SHF 


n „ i „ i 2 /Q „ i 1 — , „4 

Ci = a s + ay + ^jj s p + 2ISP 

C 2 = 5s + 5 v + 6 -fB - e -fB Q p 

C 3 = -2Bl (a s + PsP + ISP 2 ) 

C4 = — C3 —OtBq 


60 + -f P 

62 

61 


6H = -^{(VB o )-(Vp)-2(VB o )-J-6 T (VB o p)-(Vp) 
-2B 2 Vp • (2 J + ^pVp - Vp) (a s + PsP + ISP 4 ) } 



Table 2 

The coefficients of the energy density (11) for the nonrelativistic limit of RMF-PC 
compared with the corresponding terms of SHF. The RMF-PC is grouped in two 
columns. The first column collects the terms of the standard model. The second 
column adds terms stemming from tensor coupling. Terms of the RMF-PC which 
do not fit into the form (11) are collected in the last two rows. They are all related 
to tensor coupling. This table shows only isoscalar terms in either model. 



even here the form of density dependence differs. Both models yield a very 
similar description of a broad range of nuclear observables in practice. This 
hints that the differences in density dependence are all somehow compensated 
by chosing appropriate effective strengths. We urgently need observables for 
more extreme densities to pin down the differences of the models. All terms in 
the non-simple part 57i come from the relativistic tensor coupling. It is very 
hard to assess their practical importance and relative weight. One can estimate 
that they are at most of the order of the Ap terms. A detailed analysis remains 
a task for future research. Last not least, it ought to be mentioned that some 
SHF functionals carry a term oc J 2 which is not present in the above form. It 
would appear in the nonrelativistic limit of RMF-PC only in the next higher 
order. 



To summarize: The non-relativistic limit of RMF-PC recovers the basic struc- 
ture of terms in SHF. The latter is more general in that the kinetic and spin- 
orbit terms have independent parameters while they are more or less linked 
in RMF. The RMF-PC, on the other hand, adds density dependence to each 
one of the terms while SHF uses it only in the leading term. 
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4 Density dependent spin-orbit terms 

This section continues the discussions of the non-relativistic limit with a par- 
ticular emphasis on the spin-orbit potential. It has has the general structure 

W q = b 4: Vp + b' 4 yp q + c 1 J + c[j q , (12) 
where q = p,n for proton or neutron density. The last two terms are not 
found in the non-relativistic reduction. They are discarded in the following 
discussion. The parameters 6 4 and b' A are generally density dependent. The 
derivation from the RMF-PC model (with Hr taken into account) yields in 
detail 

A(po) 

{2m q + A'(p )p G + Bp q f' 
B 

(2m q + A'(p )p + Bp,) 2 ' 

A(p ) = (a s - a v + a rV ) + 2/3 s p + 3(7s - 7v)Po, 
A'(p ) = (a s -«v + a rV ) + Pspo + (7s - 7v)Po, 

B = -a rY . (13) 

A tensor term (as it appears for example in the model of Rusnak and Furn- 
stahl [22]) will induce a further correction in b 4 and b' 4 . 

The spin-orbit term which is usually employed in Skyrme energy functionals 
exhibits two significant differences compared to the nonrelativistic limit of rel- 
ativists models, namely (a) the restricted isovector dependence and (b) the 
lack of nonlinear terms. An extension of the Skyrme model with an enhance- 
ment in W by isospin contributions has already been done (SKI3-4) [17]. By 
introducing different isospin contributions into the spin-orbit potential, SKI3- 
4 do reproduce the isotope shift of the rms radii in the heavy Pb isotopes, but 
these parameter sets still yield different shell closures in superheavy nuclei 
than the RMF models [16]. On the other hand, in relativistic point-coupling 
models, the role of the density dependence has proven to be important to re- 
produce acceptable single particle spectra and spin-orbit splittings [2]. There- 
fore the enhancement of SKI3-4 with nonlinear terms (density dependence) 
might result in an improvement of their shell structure predictions. 



b 4 = 

K = 

with 
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In Fig. 1, we show the neutron spin-orbit potential for the two RMF forces PC- 
Fl (RMF-PC) and NL-Z2 (RMF-FR) as well as for the Skyrme interactions 
SkI3 and SkI4. (The spin-orbit potential for the RMF was obtained from 
mapping the Dirac equation into an effective Schrodinger equation.) Shown 
are the results for the four doubly-magic systems 16 0, 48 Ca, 132 Sn and 208 Pb. 
We recognize that the potentials of the two RMF models have a similar radial 
dependence, the potential of NL-Z2 being a bit deeper in three cases. The spin- 
orbit potential of the Skyrme forces, however, is both shifted to larger radii and 
also deeper in comparison with the RMF results. As could be shown in Ref. [2], 
the RMF-PC model with PC-LA has similarly a too deep potential peaked at 
a too large radius. It suffers from the same wrong trend with mass of spin-orbit 
splittings as do SkI3 and SkI4 (splittings get too large with increasing mass). 
The reason there are the actual values of the nonlinear parameters, leading 
to a density dependence of the mean fields that lead to this situation. This 
issue has been cured with the introduction of PC-F1 which has a predictive 
power of spin-orbit splittings comparable to the best RMF-FR forces. The 
Skyrme force SkI3 has a spin-orbit term that mimics the isospin-dependence 
of the RMF model, while SkI4 has an additional free parameter. However, 
the spin-orbit potential of SkI4 lies closer to the RMF results. A quantitative 
calculation using the Skyrme model not only with isospin terms but also with 
density dependent terms in the spin-orbit potential is strongly suggested. A 
density-dependent ansatz in the spin-orbit potential has been introduced in 
Ref. [23], but unfortunately this reference does not give a parameter set for 
this type of Skyrme model. The form of the spin-orbit potential, which is 
rather an ad hoc ansatz, is quite different from the one presented here. The 
spin-orbit part of the energy density which reproduces the spin-orbit potential 
predicted by our analysis is given by 



£is = 6 4 pV J + b' A {p p WJ p + p„V J n ) 

+\{W lP + W 2 p 2 )(J n ■ V Pp + J p • Vp n + ^2,J q - Vp q ). 

(14) 

The values of above parameters can be illustrated by taking the PC-F1 pa- 
rameter set, using p eq = p nm (PC-F1). We obtain 6 4 = 122.93 MeV fm 5 , 6^=7.01 
MeV fm 5 , W lPcq /2=3QM MeV fm 5 , and W 2 (p cq /2) 2 =-10.2Q MeV fm 5 . 
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5 Modifications due to exchange terms 



Taking into account the exchange correction in the nonlinear terms (for more 
details about the calculation of these terms see Refs. [3,19]) , after a straight 
forward calculation in the effective Hamiltonian except C 2 , all appear to be 
density dependent (since the explicit form does not provide particular insight, 
we omit presenting it here). If we calculate them by using the model of Rusnak 
and Furnstahl [22] with the derivative nonlinear terms taken into account, 
we will obtain C2 also to be density dependent. If we take into account the 
exchange corrections, we will have not only modifications in 64 and b' 4 , but C\ 
and d x also are not zero, due to the contribution of the tensor terms which come 
from the Fierz transformation [3,19], these corrections are small, however. 

The explicit treatment of exchange terms, due to Fierz transformations both 
in Iso and Dirac space, gives birth to a variety of additional isovector terms 
without introducing new parameters. This is quite interesting, since there is 
still a problem in the isovector channel of the RMF model: on one hand, it 



13 



seems to be not flexible enough to account for isovector data, on the other 
hand, modern fitting strategies fail to fix additional terms corresponding to 
extensions in the isovector channel [2]. Thus, modifications governed by the 
treatment of exchange might be a cure. 

6 Conclusions 

We have performed a nonrelativistic reduction of the relativistic point-coupling 
model and compared it to both the Skyrme-Hartree-Fock model and the non- 
relativistic limit of the linear RMF model with meson exchange. The motiva- 
tion was to gain more understanding about the interrelations of these different 
approaches, which, though looking quite different at first sight, appear to be 
very similar in the nonrelativistic limit. 

We found that there are some significant differences in the models, namely 
(a) there is a difference in the parametrization of the density dependence of 
the mean-field, (b) there is no explicit density dependence in the spin-orbit 
term in the SHF model. We have written down an ansatz for a Skyrme energy 
functional containing these extensions which should be studied numerically in 
the future to see if it can improve the spectral features of the SHF model. 

A complete treatment of exchange terms in the RMF-PC model would lead 
to a variety of additional isoscalar and isovector terms, which, in turn, would 
strongly influence the nonrelativistic limit. Relativistic Hartree-Fock calcula- 
tions open the possibility to a close investigation of nonrelativistic versus rela- 
tivistic kinematics and their relevance in effective models for nuclear structure 
calculations. 



A Appendix 

A.l Skyrme Hartree-Fock Model (SHF) 

The energy density functional of standard SHF is 
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the kinetic term is 



The Skryme part reads 



h 2 



e\ iD = —r. (A.2) 
2m 



.°sK = |p 2 -|Ep? + |p a+2 -|p a E^ 

+ h (pr - f) - b[ Y,{ Pq r q ~ f q ) ~ £pAp + |E PAP, 

q Z Z q 

- b 4 ( P V ■J + a-(Vxj)+ EtP 9 (V • Jq) + <Jq ■ (V X J q )\). 



(A.3) 



Additionally, the Coulomb energy has to be added. The densities not already 
defined previously are 

J q = ~\ E[^(9)V^(g) - {V^{q))Hi{q% (A.4) 

^ = E^!(9)^(9)» ( A - 5 ) 

i 

but these densities are zero in spherically symmetric systems due to time- 
reversal invariance. 



A.2 Walecka Model (RMF-FR) 



The Lagrangian density for RMF-FR is 

'free i /'free 

-"nucleoli ' meson """"coupl ~ r "^coupl ' 



where 



n /'free _i_ /'free _i_ nXni _i_ /"nonlin I \ p\ 

*~ ^nucleon ' ''-'meson ' ''-'coupl ' ^coupl ' v-™- -0 / 



Cdeon = ttiluP - m B )V, (A.7) 



C C e C son = ^(^^$-m^ 2 ) 

-\{\G» v G^-mlV»V») 
-\{\B» u -B^-mlR,.&) 

- ^ V, 

(A. 
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4ou P i = -9s^ - gvV^Y^P - 9rR, ■ i>?-fi> (A.9) 

ljV - ^d u R^r^r^ (A.10) 



2m,B 4m 
-eA l $?-tp-'fil>, 



(A.ll) 



and 



£SSr = -^ 2$3 - ^3$ 4 . (A.i2) 

A. 3 Nonlinear Scalar Meson Equation 

Here, we give an alternative approach [25] to arrive at Eq. (5a). In this ap- 
proach it becomes transparent that the convergence of the nonlinear expansion 
of the scalar density depends solely on g s and bi and not on the expansion value 
of the scalar density. The equation for the scalar meson can be written as 

(V - m 2 )$ = g sPs + 6 2 $ 2 + 6 3 $ 3 . (A.13) 

One way to solve this equation is by using an iteration procedure as follows: 
we start by choosing an initial $: 

$ o(ro) = -gs J D(\f - n\) Ps (n)d s ri , (A.14) 

where D(\fo — f{\) is the Greens function satisfying (V — m 2 )<E>o = 9sPs, and 
substitute this $ into 

(V - m^n = 9sPs + &2^_i + h®l-v (A.15) 
We can say that the iteration procedure terminates if <&k=&k+i- The result is 

§{r) = -gs J D(\f - f 1 \)p s {f 1 )d^r l 

+ / / f2(r,f 1 ,f 2 )p s (fl)p s (r 2 )d 3 r 1 d 3 r 2 

+ ••• (A.16) 

where 
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/ 2 (r,f!,f 2 ) = -g% J D(\r-f\)D(\f -f[\)D{^ -r 2 \)d 3 r' (A.17) 



/ 3 (r,n,f 2 ,f 3 )=A / D(\r-7\)D(\7-f 1 \)D(\7-f 2 \)D(\7-f 3 \)cPr' 

-9 3 s b 2 2 J jD{\r- r"|)L>(|r" -?\)D{\i" - fl\)D{\7 - f 2 \) 
■ D(\7 -r 3 \)d 3 r'd 3 r", (A.18) 

... is the contribution from terms with p s more than three. It is clear from 
Eq. (A. 16) that the fast convergence in the nonlinearity ( in the power of p s ) 
can only be obtained if the coupling constants g s and hi are small, basically 
independent of the point of expansion. 

Because D(\r" — f*]) is a distribution function, we can expand it into a delta 
function and its derivative: 

Dtfy) = -^5 3 (r-7) + ^-5 3 (r-f f ) + ... (A.19) 

m l s rrig 

If we choose to expand only up to the second term, insert it into $ and 
insert $ into the scalar Hamiltonian, we obtain the same result as Eq. (5a). 
This delta expansion, if integrated with density, has a connection with the 
Taylor expansion of the density. It can be easily understood from an artificial 
1 dimensional (ID) illustration as follows: in ID, the propagator can be written 
as 

= fo6(x) + h5\x) + .. + fj^(x) (A.20) 
where S (n \x)= d n /dx n 5(x). Then %(0) is 



/oo 
p(x)D(x)dx = f o p(0) - f lP '(0) + ... + (-l)V n) (0)- 
-oo 



Now if we expand p(x) around x=0 in a Taylor series as 



(A.21) 



p(x) = p(0) + p'(0)x + ... + ^P (n) (0), (A.22) 
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we obtain 

/oo poo roo 1 

D(x)dx + p'(0) / x J D(x)rfx + ... + p (n) (0) / -7X n D(x)da:. 
-oo J—oo J—oo TV. 

(A.23) 

If we compare both <3> , we have an alternate representation of D(x): 

00 {—\\ n , s r°° 1 
D( I ) = y , iiL{WM / -x in D(x')dx' 

It is clear that the choice of proper origin for the expansion in a; has an effect 
on how far we need to expand p(x) to obtain a good approximation. If we for 
example find x=0 as a good position, so that p(x) ~ p(0) + p'(0)x, it means 
that for every f n , n>l gives contribution zero. The last equation tells us that it 
is nothing else than D(x) ~ - ^5{x) +^5 f (x). We suspect a similar behavior 
to happen in the real word (3D). 



(A.24) 



A. 4 Details of the nonrelativistic reduction 



This appendix provides a brief outline of the nonrelativistic expansion of the 
scalar and tensor densities. For simplicity, the wavefunctions are used without 
the index for the state a. 

The Dirac equation (7a) decomposes into 

{m-e + S + V )^ up) + a-{p + iT)^ =0 , (A.25a) 
(m + e + S - V )ip {dw) -a-(p- iT)ip^ = . (A.25b) 

The scalar and vector densities are 



Ps = 




2 




2 


(A.26a) 


Po = 




2 

+ 




2 


(A.26b) 



The "normal" (baryon) density is the vector density po- In the following, we 
eliminate the lower component thereby carrying forth only terms up to second 
order vap/m. 



18 



The lower component can be expressed through the upper component 



ip^ = B a-(p-iT)ip 
i 

Bo 



(up) 



m + e + S -V 2m + S - V 



(A.27a) 
(A.27b) 



The upper component is not yet the nonrelativistic wavefunction because as 
such it is not normalized. The classical wavefunction is introduced through 



I=l + a-(p + iT)B%a-(p-iT) . 
This yields the desired result for the vector density 

2 



Po 



(cl) 



(A.27c) 
(A.27d) 



(A.28) 



More involved terms appear for the scalar density: 



= ^W+f-yO - <pW + i-V*a. (p + iT)B 2 a-(p - iT)t 1 / 2 ^ 
= p - 2y? (d) V (p + iT)B 2 a-(p - #> (d) 

= Po -2B 2 (t -V-J + p T 2 -2T-J + T • Vp) . (A.29) 
The tensor density is expanded as 



I —2(7 

p T = Vy? + I I if 

ia 

-V^ up)+ iaB a-(p - iT)cp {up) + V V {up)+ -a{p + iT) B iacp (up) 
-V(B V Po ) + 2V(B J)-V(B T Po ) . (A.30) 
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